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Abstract. This paper summarizes some known results about Appell polynomials and investigates 
their various analogs. The primary of these are the free Appell polynomials. In the multivariate case, 
they can be considered as natural analogs of the Appell polynomials among polynomials in non- 
commuting variables. They also fit well into the framework of free probability. For the free Appell 
polynomials, a number of combinatorial and "diagram" formulas are proven, such as the formulas for 
their linearization coefficients. An explicit formula for their generating function is obtained. These 
polynomials are also martingales for free Levy processes. For more general free Sheffer families, a 
necessary condition for pseudo-orthogonality is given. Another family investigated are the Kailath- 
Segall polynomials. These are multivariate polynomials, which share with the Appell polynomials 
nice combinatorial properties, but are always orthogonal. Their origins lie in the Fock space repre- 
sentations, or in the theory of multiple stochastic integrals. Diagram formulas are proven for these 
polynomials as well, even in the g-deformed case. 



are finite. Then there are at least three natural families of polynomials associated to such a measure. 
The most familiar ones are the orthogonal polynomials {-Pnj^o- ^^^^ ^ polynomial family (that 
is, Pn has degree n) such that 



for n ^ k. Two standard normalizations are to require the polynomials to be monic or to be 
orthonormal. By the spectral theorem for orthogonal polynomials, such (monic) polynomials satisfy 
a three-term recursion relation 



where G M, /3„ G IR+ are the Jacobi parameters, and P_i = 0, Po{x) = 1. The Jacobi parameters 
and the moments of the measure can be expressed in terms of each other, and their properties related 
to the properties of the measure and the orthogonal polynomials, for example using the Viennot- 
Flajolet theory lE^H^lHTll . A typical question in this direction is to find explicitly the linearization 
coefficients 



1. Introduction 



Let /i be a probability measure on the real line all of whose moments 





Pn+l{.x) = xPn{x) - a^Pnix) - /3„P„_i(x) 
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Already the proofs of the positivity of these coefficients are quite subtle lfT9ll and they are known 
explicitly only in very rare cases [|^ . 

Another natural and very classical fTO*! polynomial family associated to fi is its family of Appell 
polynomials, which have the exponential generating function 

oo ^ 

(1) = exp(xz - logM^(^)), 



n=0 



where 



oo ^ 



n=0 

is the exponential moment generating function of /x. It is easy to see that an equivalent definition is 
via a recursion relation 

(2) = xA^-\x) - J2 (^)r„+i_,(/i)A('=)(x) 

A:=0 ^ ^ 

and A^'^\x) = 1. Here, r^in) are the cumulants (semi-invariants) of the measure. Unlike the Jacobi 
parameters, the cumulants are obtained from the moments of the measure via a simple relation 



oo ^ 



k\ 

k=l 

(both sides considered as formal power series). Many classical (non-orthogonal) polynomial fami- 
lies are Appell. They arise in finite operator calculus iHOl and the study of hypergroups, in numerical 
analysis (Bernoulli polynomials are Appell), but also in probability theory, in the study of stochastic 
processes |32|, non-central limit theorems L13t.25 l. and natural exponential families [,39J . From the 
combinatorial point of view, they have nice linearization and multinomial formulas. 

The third family of polynomials has not apparently been explicitly defined before, although it ap- 
pears implicitly in the paper ll43l . For this reason, we will call them the Kailath-Segall polynomi- 
als. These are polynomials in (infinitely many) variables {xk}'^^i- They are indexed by all finite 
sequences of natural numbers u = (m(1), m(2), . . . , u{n)), n > 0, and defined by the recursion 

(3) W, 

,u(2),...,«(n) 
n n 

i=l 1=1 

with initial conditions = 1, Wi = Xi — r^. As usual, u{i) means "omit the i'th term". 
Note that Wn contains a single monomial of the highest degree \u\ = ^"=1^(0' namely = 
Xu{i)Xu{2) ■ ■ ■ Xu{n), and that it is a polynomial in the variables 



u{j) for some S C {1,2, ... ,n} > 
ies J 



Xi 

The origin of these polynomials is again in probability theory, where they appear as certain multiple 
stochastic integrals. Moreover, they share a number of properties with both the orthogonal and 
the Appell polynomials. Their recursion relation is determined by the cumulants, and they have 
nice linearization properties, just like the Appell polynomials. In fact, we will show that A^'^\xi) 
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is a linear combination of the Kailath-Segall polynomials. On the other hand, the polynomials 
= 1 1 also have the following orthogonality property. Define a measure /x*^") on M" by 
specifying its multivariate cumulants: 

See Remark [U for the fashion in which the measure /i*^"-' is determined by its cumulants. If fi is 
infinitely divisible, fi^'"^ is a positive measure. Then 

/ W^"'\xi,X2,...,Xn)W^''\xi,X2,...,Xk)dfi^''\xi,X2,...,Xn) = 

if k < n. 

A generalization of the Appell polynomials are the Sheffer polynomials. Let f/ be a function with 
a formal power series expansion such that U{z) = z + higher-order terms. Then the Sheffer 
polynomials are defined via their exponential generating function 

oo ^ 

-^Pn{x)z'' = exp(xf/(^) - \ogM^{U{z))). 

n=0 

Equivalently, 

U-\d,){Pn) = nP^_i and (P„, 1)^ = forn > 0. 
They share with the Appell polynomials the multinomial expansion properties and the relation to 
stochastic processes (see Section l2!4l) . Among the Appell polynomials, only the Hermite ones are 
orthogonal. Meixner's classic characterization 0?ll describes all the orthogonal Sheffer polynomi- 
als. There are also multivariate versions of this statement [38|. 

We start the paper by describing some properties of and relations between the three aforementioned 
families of polynomials. In the Appell and Kailath-Segall cases, natural starting points are in fact 
certain multi-linear functionals on more general algebras, which can then be specialized to polyno- 
mials. These families have similar but different properties. For example, the Appell linearization 
coefficients are sums over non-homogeneous partitions, while the Kailath-Segall linearization co- 
efficients are sums over inhomogeneous partitions. A number of the results in SectionElare known 
and so are stated without proof. 

As mentioned above and described in more detail below, both the Appell and the Kailath-Segall 
polynomials arise in probability theory and are related to the notion of independence. Let be a real 
linear functional on the algebra Mfx, of polynomials in two variables. Then x,y are independent 
with respect to if for any P, Q, 

if[Pix)Q{y)]=cp[P{x)]ip[Qiy)]. 

In the early 1980's, Dan Voiculescu introduced a parallel but really very different notion of free 
independence ll48l . Let tp now be a real linear functional on the algebra M(x, of polynomials in 
two non-commuting variables. Then x, y are freely independent with respect to Lp if whenever 

[Pl(x)] = ip [QM] = ... = ^ [Pn{x)] = ^ [Qn{y)] = 

and Qo, Pn+i each are either centered or scalar, then 

if [Qo{y)Pi{x)QM ■ ■ ■ Pn{x)Qn{y)Pn+i{x)] = 0. 

A whole theory of free probability based on this notion, is by now quite well developed. It 
turns out that there are "free analogs of" the Appell and Kailath-Segall polynomials, which, very 
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roughly, are obtained by replacing commuting variables with non-commuting ones, exponential 
generating functions with the usual ones, and exponentials with resolvents. Such a replacement, 
however, is quite non-trivial. The analysis of the preceding section is repeated in Section|5]for the 
free case, except that in this case most of the results are new. For the free Appell polynomials, 
we find an explicit form of the generating function and various recursion and "diagram" formulas. 
As expected, these formulas are in most cases based on the lattice of non-crossing, rather than all, 
set partitions (but not always, compare for example equation (fT3t with Proposition fT2l c)). For 
the free Sheffer polynomials, we find a necessary condition for them to be pseudo-orthogonal. 
Consequences of these results for free probability will be developed elsewhere. 

Comparison of formulas from the preceding two sections shows that many of them appear as par- 
ticular cases of g-interpolated forms, with the usual case corresponding to g = 1 and the free case 
corresponding to g = 0. On the other hand, many other formulas do not appear to admit of such an 
interpolation. In Section|4l we define the g- Appell and g-Kailath-Segall polynomials, and show that 
some of their properties carry over to the whole interpolated family. Some of the other properties, 
at least at present, do not. Therefore this section is necessarily more tentative than the preceding 
ones. We only consider single-variable g- Appell polynomials, and find an explicit form of the gen- 
erating function for them, as well as the g-analogs of the Meixner families. For the g-Kailath-Segall 
polynomials, we find various recursion and "diagram" formulas. Finally, in the appendix we show 
that g-Appell polynomials are not linear combinations of the g-Kailath-Segall polynomials, unlike 
in the classical and the free case. As a consequence, they cannot be martingale polynomials for the 
g-Levy processes. 

Acknowledgements: Thanks to Mourad Ismail, Marius Junge, and Michael Skeide for useful and 
enjoyable conversations. Thanks also to the referees for a number of helpful suggestions and criti- 
cisms, especially the correction in Definition |5] 

2. Classical polynomial families 
2.1. Notation. We will use multi-index notation u = (m(1), m(2), . . . , u{k)). Denote by 

^k,n = {u e N''' : \u\ = n] 

a basic simplex (where ^ N). For two multi-indices u, v, {u, v) will denote their concatenation. 

For u G N", 5 C {1, 2, . . . , n}, 5 = {v{l) < v{2) < ... < v{k)}, denote 

{u:B) = {u{v{l)),u{v{2)),...,u{v{k))). 
For a subset A C B, will denote the complement of A, where B is understood. 

2.1.1. Partitions. A set partition of a set S' is a collection of disjoint non-empty subsets of 5* whose 
union is S. If S is an ordered set and tt = i?2, • • • , -Bfc} is such a partition, we order the classes 
of TT according to the order of their smallest elements, min(i?i) < min(i?2) < . . . < min(i?fc). 
We will consider the following three lattices of partitions. By V{n) we'll denote the lattice of 
all partitions of the set {1,2,..., n}. For tt G V{n), Sing(7r) is the collection of single-element 
(singleton) classes of tt. 

By NC{n) we'll denote the lattice of non-crossing partitions iBTl . These are the partitions with the 
property that 

TT 

i<j<k, i'^k, j'^l, ir^j =^ i < I < k. 
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For a non-crossing partition vr, a class B is called outer if for any other class C G vr, if i,j G C, 
k E B, and i < k, then j < k. Otherwise a class is called inner. The outer classes of tt will be 
denoted by Outer (tt). 

The third lattice of partitions, which we use mostly for notational convenience, is that of interval 
partitions, all of whose classes are intervals of consecutive integers. This lattice Int{n) is naturally 

isomorphic tothe power setof {1,2, ... ,n — 1}. For^^^j^s(i) = n, we denote by 7rs(i)^s(2),...,s(fc) ^ 
Int{n) the interval partition 

(l,...,s(l)),(s(l) + l,...,s(l) + s(2)),...,(5^s« + l,...,n) I. 

i=i J 

There is a partial order < on V{n) which restricts to the other two lattices. We denote the smallest 
element in V{n) by = {(1), (2), . . . , (n)} and the largest one by 1 = {(1, 2, . . . , n)}. We denote 
the meet and the join in the lattices by A and V, respectively. In particular, 

I ~ J <^ 1 ^ J and « ~ J. 

Definition 1. Let a G V{N) be a partition. For a partition tt G V{N), we say that 

(a) A class i? G vr is homogeneous with respect to a if i? C C for some C G a, 

(b) TT is non-homogeneous with respect to a if vr has no homogeneous classes with respect to cr, 

(c) vr is inhomogeneous with respect to a if tt A a = 0, that i^, i j ^ i j , 

(d) vr is connected with respect to cr if tt V a = 1. 
See Figure [2 



a) 




Figure 1. With respect to the partition {(1, 2, 3), (4, 5, 6), (7, 8, 9)}, partitions 
which are (a) non-homogeneous, connected, (b) inhomogeneous, connected, (c) not 
connected. 

2.1.2. Extended partitions and restricted crossings. For S <Z -k, call the pair (S*, vr) an extended 
partition; S is to be thought of as the collection of classes of tt "open on the left". See FigureElfor 
an example. 

For 1 < k < m < n, define the restriction 

{S'y) = {S,Tx) Kfc,...,m} 
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Figure 2. An extended partition of 10 elements with 2 left-open classes and 4 
restricted crossings. 

as follows: 

B' e 7t' if B' = B n {k, . . . ,m} , B e 7T, 
B' e S' if B e S or B n {1, . . . ,k - 1} ^ dl. 
Define the number of right restricted crossings of (S, n) at the point k as follows: 

{0, if keB,k = max{B) ork+leB, 
\S'\, if k e B,j = min{i e B,i > k} , 
iS',n') = {S,n) f {A; + l,...,j-l}. 

Let rc {S, vr) = X]fe=i {k, S, vr). Note that also 

10(5, tt) = rc (tt) + ^ |C G tt : min(C) < min(i?) < max(C)| , 
Bes 

where rc (vr) = rc (0, tt) (see [14]). 

2.1.3. Cumulants. A measure on M all of whose moments are finite induces a positive semi- 
definite unital linear functional Lp on Wyx] by Lp [x"] = mn{iJt)- Positivity will not play a part in most 
of the results below. Thus, our starting object is a linear functional Lp, which does not necessarily 
correspond to a positive measure, although we still assume that it is unital, Lp\\] = 1. Throughout 
the paper, the functional will be fixed, and so will frequently be omitted from notation. 

In the multi-dimensional situation, let be a general unital real linear functional on some real 
algebra A. For Xi, X2, . . . , X„ G ^, denote by 

M[Xi,X2,...,X„] = (/.[XiX2...X„] 

the joint moment of (Xi, X2, . . . , X„) with respect to ip. Also, for tt G V{n), denote 

M^[Xi,X2, . . . ,X„] = J] M[X, -.zeB] 

Ben 

the partitioned moment of (Xi, X2, . . . , X„). Here, in the case when A is non-commutative, the 
factors in Ylies terms in {Xj : i G B} are taken in order. For each vr, is an n-linear 

map. Given M, define the corresponding jom? cumulant recursively by 

(4) _R[Xi, X2, . . . , X„] = M[Xi, X2, . . . , X„] — ^ _R^[Xi, X2, . . . , X„], 

■KeV(n) 



APPELL POLYNOMIALS AND THEIR RELATIVES 



7 



where the partitioned cumulants Rt, axe defined as above, and involve only cumulants of lower 
order. These are also multi-linear maps. 

The multi-index notation will be used extensively but consistently. Throughout the paper 



On the other hand, x = {xi, X2, ■ ■ ■ , Xn) will denote the total collection of variables involved, and 
same for z, w, X, etc. Finally, for i? C {1, . . . , n}. 



the product taken in order of increasing indices. 

For S = {Bi, . . . , Bk] a collection of disjoint subsets of {1, ... , n} (for example, a partition), S is 
always ordered according to the order of the smallest elements of the subsets. 

Remark 1. If is a measure on R" all of whose moments 



are finite, its cumulants r^(yu) are defined in terms of the moments by equation ©. Note that while 
the functional xu ^— ^ rn^^jj) is determined by its cumulants, again via equation Q, the measure /z 
itself may not be. The determinacy of the moment problem in the multivariate context is a difficult 
question, see [20, Chapter 3] and their references. 

2.2. Appell polynomials. Appell polynomials are defined by equations O or Q. They satisfy the 
following properties, easily obtained using the generating function: 



(5) 




XB = \{ 




(6) 



(7) 




k=0 



n 



(8) 




k=0 



where 




and the notation is suggested by the relation 




See the original paper ifTDIl . where ® was taken as the definition, or lE51l . 
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Remark 2. The following result of Appell ifTOll . developed in great detail by Rota et al. f^O"], appears 
to have no free analog. Let { A^") } , | } be two families of Appell polynomials, with exponential 
generating functions F, G, respectively. Define a new family of polynomials { (AB)^"^^ } as follows: 
expand A*^"^ in the powers of X, and for each substitute B^^l Then = and 

these polynomials are again an Appell family, with exponential generating function FG. 

2.3. Multivariate Appell polynomials H^H^Hl. Let ^ be a commutative real algebra with a unital 
real linear functional if. We will usually call elements of such an algebra random variables, since 
any collection of real-valued random variables on some probability space, such that all of their 
joint moments are finite, generates a commutative real algebra with the expectation functional on 
it. For n > 1, define an n-linear map A : — > ^ as follows. For Xi,X2, . . . G A, 
A (Xi, X2, . . . , Xn) is in fact a symmetric polynomial in {Xj}"^^, which we denote by 

It is determined recursively by 

(9) c?x,^Xi,X2,...,x„ (a^i, a;2, . . . , x„) = ^Xi,X2,...,x,,...,x„ (^i' • • • , ^i, • • • , a^n) 
and 

(10) <^[A(Xi,X2,...,X0] = 

for n > 0, where = 1 is understood. The advantage of this notation is that we can consider only 
the algebra generated by {X^}"^^ (and 1), and the restriction of ip to this algebra can be thought of 
as the joint distribution of these random variables. For u E {1,2,..., nY, define the polynomials 

Au (x) by 

An (x) = ^X„(i),X„(2),...,X„(fe) X„(2), . . . , Xu{k)) ■ 

Thus 

(X) = A (Xu(i), Xu(2), . . . , X„(fc)) . 

Note that this notation differs from the usual one: A^ depends in fact only on the number of occur- 
rences of each index u{i), and not on their order; in the usual notation one writes down the number 
of such occurrences. Our notation is better suited to the non-commutative case. 

The following are some properties of the multivariate Appell polynomials. Denote 

n 

X ■ Z = ^ ^ X^Zi, 
i=l 

and denote by 

00 _ 

^==1 v€{l,...,n}'' 

the exponential cumulant generating function. Note that the 2j's commute and R is symmetric in its 
arguments, so each term 

n 

appears 

z{2), t{n)) ~ ^(l)H(2)!...^(n)! 
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times. So this notation coincides with the usual one. 
Single- and multivariate Appell polynomials are related via 

v4(")(x) = Ai,i,...,i(a;,x,...,x), 
where the single-variable polynomial corresponds to the moment sequence m„ = y9 [X"]. Also, 



(11) 



l + ^ry X] A^(x) 2;^ = exp(x ■ z - R(z)), 



^==1 »/G{l,...,n}'= 



(12) A{X„X^)=X,AiX^)- R[X„X^^.,v)]A{X^^.,vc)), 

Vc{l,...,n} 

TTeV(n) Be-TT CGtt, Vc{l,.--n} TrG'PCV':) 



(13) A{x^)= j2 E^(-s)(-i)'"'"' n^[^(- 



C) 



Here is the complement of V in {1, . . . , n}. Equation (fTTT) looks slightly unusual because of 
our different notation. The single- variable polynomials also satisfy a binomial formula: if X, Y are 
independent, 



(14) {X + Y)=J2 (^) ^^^'^ (X) 

i,_n \ / 



k=0 

This is, of course, the short-hand notation for 



1,— n V / 



\x)A'r'\y)- 



k=0 

Remark 3. For fixed {Xi,X2, . . . ,Xn} C A, their joint distribution is the functional (/9x on 
R{xi,X2, . . . , Xn) determined by 

Any such functional ip on M(xi, X2, ■ ■ ■ , x^) can be included in a one-parameter family {V^tj^gjQ q^-, 
in the following fashion. Let R denote the joint cumulants of {si, X2, ■ ■ ■ , Xn} with respect to if). 
Then define 

i^Axu] = E n R{^{u--B)]- 

That is, the joint moments of {xi,X2, ■ ■ ■ , x„} under ipt correspond to the joint cumulants tR via re- 
lation Q. If all the linear functionals ipt are positive, the corresponding measures form a semigroup 
with respect to the convolution operation. They are also marginal distributions of the corresponding 
Levy process (see Section l2!4b . All {tpt} are positive if and only if tp itself is infinitely divisible. 
Without the positivity requirement, any functional is algebraically infinitely divisible. So any family 
of Appell polynomials naturally comes included in a one-parameter family A^'^\-, t). 

A similar construction, based on relation (fT9b can be done in the free case (see Section 13). In this 
case, if all the linear functionals 'ipt are positive, the corresponding measures form a semigroup 
with respect to the additive free convolution operation, and they are marginal distributions of the 
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corresponding free Levy process. This is the case if ip itself is freely infinitely divisible. However, 
in the free one-dimensional case, more is true: for any positive tp, is positive for t > 1 [35 1. 

2.4. Martingales. Let {X{t)}^^^Q be a Levy process, that is, a stochastic process with stationary 
independent increments. Assume that all of the joint moments of {X(t)} (with respect to the 
expectation functional) are finite. Let be the conditional expectation onto the (von Neumann) 
algebra generated by {X(s)}^g[p which extends to an orthogonal projection on the space of all 

square-integrable random variables. Then for each n, A^'^\X{t)) is a martingale, that is, 

E, [A(")(X(t))] =A(")(X(s)). 

The following are two elementary proofs of this fact. If M{z) is the moment generating function for 
X(l), then the moment generating function for X(t) is M^zf. So using the generating function (HI) 
of the Appell polynomials, 

E, [exp(X(t);z - tlogM(^))] 

= E, exp(^(X(t) - X{s))z -{t-s) logM(z)) exp(x{s)z - shgM{z)^ 
= exp{X{s)z - s\ogM{z)) 
since E [exp((X(t) — X(s))z)~\ = M{zY^^ . On the other hand, using the binomial property (fT4l) . 

E, [A(")(X(t))] = E, U(")(X(.) + {X{t)-X{s)) 



E. 



k=0 ^ ^ 



A(")(X(.)), 



since E [A*^" ^\X{t) — X(s))] = forn > A;. Independence of increments and properties of the 
conditional expectation are used in both proofs. 

Definition 2. A monic polynomial family in n variables is a family of polynomials indexed by 

oo 

ue[_}{l,2,...,nf 

k=0 



such that P0(x) = 1 and each 

Pu{x) = Xii+ lower order terms, 

where the grading is by total degree in xi, X2, . . . , x„. All the families of polynomials considered 
in this paper form monic polynomial families. 

Let M(x) = M(a;i, X2, . . . , Xn) be all the real polynomials in n non-commuting variables. For a 
multi-index m G N^, denote 

iur = iuik),...,ui2),uil)). 
Define an involution on ]R(x) via an M-linear extension of 

Similarly, define an involution on C(x) via a C-anti-linear extension of the same relation. 
A monic polynomial family {Pu} is pseudo-orthogonal with respect to a functional ip if 

^ [P*uPv] = 
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whenever \u\ ^ \v\. The family is orthogonal if this is the case whenever u ^ v. 

2.5. Fock spaces. Let {P^(X)} be a monic polynomial family in n variables. Define a func- 
tional if on C(Xi, X2, . . . , Xn) by (p [1] = 1, (p [Pn] = for \u\ > 1, and extend C-linearly. On 

,X2,...,Xn), define an inner product via 

{Xit, x^) = (f [P^Pi ■ 

So this is nothing other than the Gelfand-Naimark-Segal construction. Note that this inner product 
need not be positive; it will be positive (resp, positive definite) if the functional (p is. Define the 
action of C(Xi, X2, . . . , X„) on C(xi, 0:2, . . . , x„) by 

Pu ' 1 ■^u 

and 

fit ' i^PuPv) ' !• 

Then for any P, 

y.[P] = (l,P-l). 

Pu(X) is sometimes called the Wick product of X„(2), . . . , X^^k), and denoted : X^ :. 

2.6. Appell Fock space. Given a family of polynomials (xi, 0:2, . . . , Xn) satisfying 

j:u(j)=i 

there is a unique functional if on M[xi, X2, . . . , x„] such that Lp [1] = 1, ip [A^i (x)] = for all u. 
The Fock space construction provides a natural way to recover such a functional. The induced inner 
product on C(xi, 0:2, . . . , x„) is determined by 

l/C{l,.--,A:},C/^0 
Vc{l,...,l},Vj^9 

for u E N^, {7 G N'. Since in this section the cumulant maps R are symmetric in their arguments, 
the inner product is degenerate and factors through to C[xi, 0:2, • . . , Xn\- 

From the recursion relation (fT2t . the action of the operator Xj is 

XjXit = XjXit+ P[Xj, X(^.y)]X(t/;yc). 

Vc{l,...,n} 

Thus it is a sum of a creation operator 
a scalar operator 

and some unusual annihilation operators. 
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2.7. Sheffer and Meixner families. Sheffer families are monic polynomial families {P„ (x, t) 
such that their exponential generating function has the form 

oo 

H{x, t,z) = y^ -Pn{x, t)z- = /(z)*e^(^)^ 

n=0 

for U{z) = z + higher-order terms, f(z) = 1 + higher-order terms. Among these, the Meixner 
families are those consisting of orthogonal polynomials (see ll42l Chapter 4] for references). It 
follows from the results of Meixner |34| that up to affine transformations, the polynomials from 
this class satisfy recursion relations 

xPn{x, t) = Pn+i{x, t) + anPn{x, t) + n(t + h{n — l))P„_i(x, t) 

for a E M., h E Here we have assumed, without loss of generality, that Pi(x,t) = x. By 
re-normalizing, we can restrict the values of the parameters to the following five classes, labelled 
by the names of the corresponding families. 

Hermite: a = b = 0. Orthogonal with respect to the Gaussian distribution. 
Charlier: b = 0, a = 1. Orthogonal with respect to the centered Poisson distribution. 
Meixner-Pollaczek: 0<a<2, 6=1. Orthogonal with respect to the centered continuous 
binomial distribution. 

Laguerre: a = 2,b = 1. Orthogonal with respect to the centered Gamma distribution. 
Meixner: a > 2, b = 1. Orthogonal with respect to the centered negative binomial distribu- 
tion. 

Here for any measure fi with mi(/i) < oo, the corresponding centered measure is 

HciS) = fi{S + mi(/i)), 

for which 

xdfic{x) = mi(/ic) = 0. 

2.8. Kailath-Segall polynomials. Basic Kailath-Segall polynomials have appeared in the 
paper [43 1. They were defined via certain stochastic integrals, but the authors also showed that they 
satisfy a recursion 

(15) = V(-l) \ ^",,. a:fc+iW^^"~^^ - nW^''^ 

[n — ky. 

(with a different notation and slightly different normalization). We define multivariate Kailath- 
Segall polynomials by recursion Q; then = Note that this is a natural definition 
in the compound Poisson ("de Finetti") case; there is an analogous definition for the more gen- 
eral ("Kolmogorov") case of functionals with finite variance, see 0. That paper also details the 
stochastic integral connection in Section 3.1. 

The following more general definition comes naturally from a Fock space construction, see Sec- 
tion |431 

Definition 3. Let be a commutative complex star- algebra without identity, and (■) a star-linear 
functional on it. Denote by ^q'^ the self-adjoint elements of the algebra. Let A be the complex unital 
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Star-algebra generated by commuting symbols {X(f) : / G ^o''} (^^^ 1) subject to the linearity 
relations 

X{af + (3g) = aX{f)+f3X{g). 

Equivalently, A is the symmetric tensor algebra of ^o- The star-operation on it is determined by 
the requirement that all X{f), f G A^°' are self-adjoint. For such fi, define the Kailath-Segall 
polynomials by W (/) = - (/) and 

n 

W(/,/i,/2,...,/n) =X(/)W (A, /2,...,/n)- 5^ (//^)W (/!,...,/„...,/„) 

i=l 

n 

- 5^ W (//„...,/„...,/„) -(/)W(/i,/2,...,/n). 

1=1 

The definition can be extended in a C-linear way, so each W is really a multi-linear map from 
to A, which turns out to be symmetric in its arguments. 

In the particular case = Co[x] (polynomials without constant term), we may denote Xi = X(x*). 
The functional can be taken to be the cumulant functional of a measure /i, (x*) = ri(n). Then 

are multivariate polynomials in {xj : i G N}. If is the functional on C[xi, X2, . . .] determined by 

<^ [W^{^)] = 
for all non-empty u, then its cumulants are 

rkixn) = r„fc(/u), 

and more generally 

R[xu] = r\a\{fi). 

If /i is infinitely divisible, the functional (■) is positive. It follows that the functional if is also 
positive, see Section |431 

Proposition 1. The following expansions hold. 

(a) Of usual products in terms of the Kailath-Segall polynomials: 

X(/0X(/2)...X(/„)= Yl E n </^)-W(/c:Cg5). 

(b) Of the Kailath-Segall polynomials: 

w (A, /2, , . . . , A) = 5^ n n (ic*! - mifc). 

■K£V(n) 5cSing(7r) BeS C&S'' 

(c) Of products of Kailath-Segall polynomials: for u = . . . , Uk), 

k 

n W (/n,(l), /«,(2), • • • , /«,(s(i))) 
i=l 

= E E n (/(«^^^)>W(/(«^:C):CG^). 

■k£V(N) SCtt, B&S" 

^A^.(i),.(2),....s(fc)=oSmg(5-)=0 
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Proof. We only consider part (c), see Theorem 1261 for the remaining proofs. We use induction on 
the length of the longest Ui. If all of them have length 1, the desired statement is 

k 

nw(/.)= E E X{UB)W{fc:CeS). 

i=l TT€V{k) SCtt, BeS'= 

Smg(5'=)=0 

This is just the sum in part (a) except that is not allowed to contain any singletons, and the result 
follows from part (a) and the fact that W (/j) = X(fi) — {fi). Now let Ui be the longest multi-index. 
Suppose ^1(1) = j, and denote v = (uj(2), Uj(3), . . .). Using the defining recursion relation, 

w (A^j w . . . w . . . w (/.^j 

= W (/^J . . . lx{fj)W iU) - J2 (/i/«.W> W (^(1), . . . , A,(Z), . . . , fuM^))^ 
^ 1=2 
n N 

{f^f^m^ • • • ' A.(0' • • • ' fu^isii))) - ifj) w iU) j . . . w (/^J 
1=2 ' 

= W(^J...W(/,)W(^)...W(^J 

s(i) 

- W (/^/J . . . ^ {fjfu,{l)) W /„,;(!), . . . , fuiil), fuMi))) • • • W (/«J 

1=2 

n 

- W (/sj . . . ^ W (^fjfu.H), ■ ■ ■ , fu,{l), fuMi))) ■ ■ • W (/"J 

1=2 

Apply the induction hypothesis to the right-hand-side. The desired sum on the left-hand- side is over 
pairs (S, vr), tt inhomogeneous. Any such term appears in the sum corresponding to the first term 
on the right-hand- side. We need to show that all the other elements in the sum corresponding to this 
term cancel out. Indeed, take any pair S C n,n inhomogeneous with respect to 7rs(i),...,i,s{i)-i,...,s(k), 
but not with respect to 7Ts{i),...,s{i),...,s{k)- Any such partition contains a class B which contains the 
position of and a position of some other Ui(l). If B ^ S and \B\ =2, it gets cancelled by the 
corresponding term from the second term on the right-hand- side. Otherwise it gets cancelled by the 
third term on the right-hand-side. □ 

Corollary 2. The linearization coefficients for the Kailath-Segall polynomials are sums over inho- 
mogeneous partitions with no singletons: 

k 



n W (/„,(!), fu,{2), • • • , fu,{s{i))) 



i=l 



E n (/<»->) ■ 



TTeV(N) B£n 

tAtTsJi) 3(2) s(fe)=0, 

Sing(7r)=0 



Corollary 3. are pseudo-orthogonal. 

Proof. For m G N", {7 G N'^, by the preceding corollary 



(TeSym(n) i=l 
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Now assume that (■) is positive, (x*) = rrii^i') for some positive measure v. (In this case /i 
is infinitely divisible, in fact a compound Poisson measure, and u is the Levy measure for the 
convolution semigroup {yUt}-) Let {pi}^^ be the orthogonal polynomials with respect to u, and 
Ui = X{pi). Equivalently, {yi} is the orthogonalization of {xj} with respect to the inner product 

{xi,Xj) = ri+j{fi) = mi+j{u). 

Corollary 4. By a linear change of variables, define polynomials 

Wi{^) = W^{y). 

Then are orthogonal, that is. 

Free and even g-analogs of CoroUarylUhold, derived from appropriate modifications of Corollary|2l 
These properties are related to the "generalized chaos decomposition property" of lITTll . 

Appell polynomials are linear combinations of the Kailath-Segall polynomials of the same degree. 
Note that a priori, such a linear combination is a multivariate polynomial, but in this case it turns 
out to depend only on a single variable. 

Proposition 5. For x = Xi, 



x] 



Tr£-P(n) 
■K={Bi,B2,...,Bk} 



k=l SeA 



k,rL 



where = |{vr G V{n), tc = {Bi, B2, . . . , Bk} , u{i) = \Bi\ ,i = 1,2, . . . , k}\. Since in this case 
the variables commute and is symmetric in its indices, also 



' ^ 1 



(l!)pi . . . ...pj 



W^pi l's,p2 2's,...,pn n's(x), 



where p\- n is a number partition of the number n, ^11=1 Wi = ^^ 

Proof. The Fock representation of Section IT51 is clearly faithful. So it will suffice to show that in 
the Fock representation of the Kailath-Segall polynomials. 



• 1 



E 



Bi|,|B2|,.-,|Bfcl 



fx)-l 



7r={Bi,B2,.-,Bfe} 



Y n^i« i 



The proof will proceed by induction. Using the recursion relation for the Appell polynomials and 
induction. 



(16) 



A(-^)(x).i= Y x-U^iB.i-Yiiy^^^-^ Y U^iB 

■KeV(n) i=l k=0 ^ ^ aeV{k) i=l 
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The action of the operator X is determined by the recursion relation Q for the Kailath-Segall 
polynomials, 

X-x^= (XiW^(X)) ■ 1 = Wi^u{i)M2\...Mk) 



i=l 



+ ^ ] ^l+u{;i)Xu{l) ■ ■ ■ Xu{i) ■ ■ ■ Xu{k) + ^ ] . . . Xu(i) ■ ■ ■ Xu{k) + 



rixu- 



1=1 



Therefore expression ([T6b is 



E 

■neP(n) 





TT 


Xi ^ 






= 1 



^lX\B,\...X\B,\...X\ 



i=l 



^|B,|+l3;|Bi| • • • X\Bi\ ■ 



Xi 



+ ri 



^B, 



4 = 1 



1 = 1 



fc=0 



The first term in the sum, as well as the fourth term, are indexed by all partitions of (n + 1) whose 
first class is a singleton. The second term, as well as the third term, are indexed by all partitions of 
[n + 1) whose first class is not a singleton. Finally, the term which is subtracted is indexed by all 
partitions of + 1), with the cumulant factor corresponding to the first class of such a partition, and 
the binomial factor accounting for the choice of the remaining n — k elements of that class (other 
than the element 1). Therefore the third and fourth terms in the brackets cancel exactly all the terms 
that are subtracted, and the first and second terms add up to the desired sum over V{n + 1). □ 

Example 4. Hermite polynomials are orthogonal, Appell, and Kailath-Segall (with Xi = for 
i > 1). Charlier polynomials are orthogonal, Sheffer, and Kailath-Segall (with all Xi = x). Jacobi 
polynomials are orthogonal, Bernoulli polynomials are Appell, Abel polynomials are Sheffer. 
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3. The free analogs 

3.1. Notation. 

3.1.1. Non-commutative power series. In this section we will use power series in non-commuting 
variables. Most theorems about formal power series remain valid in this context. In particular, 
series with non-zero constant term have inverses with respect to multiplication: 

(1 - ^ auzu){l + ^ h^z^) = 1 

■3 

for 

n 

(17) = ^ aabif = ^ ^ aa^aa^ ■ ■ .a^t^, 

(u,v)=w, k=l [ui,U2,...,Uf^)=w 

where w G N". From general theory, left inverse and right inverse are equal. 

3.1.2. Free cumulants. Let G{z) = X^^o"^"-^ ^ formal Laurent series, a generating 
function for a moment sequence. Define the corresponding free cumulant generating function R(;z) 
via the functional relation 

(18) g('-±^]^z; 1±^^^.. 



Giz) 

Note that we use the boldface notation to distinguish R from the usual i?-transform R(z) = R(-2) / z 
ll49ll . Define the corresponding free cumulant sequence via its generating function 

oo 

'rnZ-\ 

n=l 

Let ^ be a possibly non-commutative real algebra, and a unital real linear functional on it. For 
X = (Xi, . . . , Xn) C Aa collection of non-commutative random variables, define their joint 
moments 

M[X^] = ^ [X^] 

as before, but now 

M{z) = J2M[X^]za 

it 

denotes the ordinary moment-generating function. Here, and in the sequel, Zi's are formal non- 
commuting indeterminates. Define the joint free cumulants of X via 

(19) R[Xi = M[X^]- n^[^(-^)]' 

neNCin), BGtt 

which expresses R[Xu{i), Xu{2), • • • , Xu{n)] in terms of the joint moments and sums of products of 
lower-order cumulants. Define the free cumulant generating function via 

R{z) = Y,R[Xu]zu. 

u 

For a single random variable X, R[X, X, . . . , X] = r„ for the moment sequence [irik = ^ [X'^j }. 
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The relation between joint moments and joint free cumulants can be summarized in a relation be- 
tween their generating functions, as follows. The following proposition is due to Nica and Speicher; 
for completeness, we provide a direct proof. 

Proposition 6. Let Zi = Wi{l + M(w)). Then 

R(z) = R(^wi(1 + M(w)),...,w7„(1 + M(w))) = M(w). 

Proof. In the defining relation dT^ . the sum is over all non-crossing partitions. Each non-crossing 
partition can be described by a subset V containing 1 and a collection of non-crossing partitions on 
the intervals into which V subdivides the set {1, 2, . . . , A;}. Applying the formula dT^ again to each 
of those intervals, we obtain 

oo s 
s=l V'C{l,2,...,fc} j=l 

y={i=i(i),i(2),...,»(s)} 

whereby convention i(s+l) = A;+l and M[X^:{j(j)+i_..._j(j+i)_i}] = M[0] = lifi(j + l) = + 
So 

oo s 

u s=l v£{l,...,ny v-i,V2,...,v^; j=l 

possibly empty 

oo s 

= Y1 E ^[^-^ n Ko) (1 + ^(^))) = (1 + M(w)) ,...,Wn{l + M(w))) . 

s=l vG{l,...,n}'' j=l 

□ 

Example 5. The key distribution in free probability, which appears for example in the free central 
limit theorem, is the semicircular distribution. It is characterized by the property that R(z) = z"^, 
its moments are the Catalan numbers, and the corresponding measure has density 



— V4 - a;2l[_2,2] (x) dx. 

More generally, we say that X has a (scaled, shifted) semicircular distribution if R(z) = az + bz'^. 
A family {Xi, X2, . . . , X„} form a semicircular family if all of them are self-adjoint and all their 
linear combinations have (scaled, shifted) semicircular distributions. Equivalently, they are self- 
adjoint and their free cumulant generating function is quadratic. They form a free semicircular 
family if in addition they are freely independent, in which case it suffices to require that each Xj be 
self-adjoint and have a semicircular distribution. 

3.2. Free Appell polynomials. 

Definition 4. Free Appell polynomials are defined via their ordinary generating function 

00 ^ 

^ ^ ' l-xz + R(z) 

n=0 ^ ' 
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Proposition 7. Some properties of the free Appell polynomials are: 

n 

(20) xA^^\x) = + rn+i-kA^'\x), 

k=0 

k=0 ^ueAfc+i,„+i i=0 

k=0 ^ S=l U&A^^n-k 

See Propo sition fT2l for a more general statement. 

Example 6. Chebyshev polynomials of the second kind, Un{2 cos 6) = ^^^^^j^^j^, are free Appell. 
They are orthogonal with respect to the semicircular distribution. In fact, from equation (l20b it 
follows that these are the only orthogonal free Appell polynomials. 

3.3. Multivariate free Appell polynomials. Define the partial derivative di on polynomials in 
non-commuting variables xi, X2, ■ ■ ■ , Xnhy a linear extension of its action on monomials 

diXu ^ ^ 3^n{l) • • • ■^u{j) • • • •^u(k)- 

j:u{j)=i 

Definition 5. Let ^ be a non-commutative real algebra with a unital real linear functional ip. We 
will continue to call its elements (non-commuting) random variables. For Xi, X2, . . . , Xn G A, 
define the multivariate free Appell polynomial 

^Xi,X2,...,x„(a;i,X2,...,x„) 

by the following conditions: 

(21) d,A {Xi, . . . ,Xi, . . . ,Xn) = v4xi,...,X,_i {Xl,. . . , Xi^i) ■ [Xi+l, ■ ■ ■,Xn) 

(with Adi = 1), 

(22) if [Ax„x„...,x^ (Xi, X2, . . . , X„)] = 0, 
and the restriction that 

(23) in each monomial, the variables appear in the increasing order of indices. 
The polynomials are uniquely determined, and since for i < j, 

didjAx^,X2,...,Xn {Xi, . . . ,Xi, . . . ,Xn) 

= {Xi,. . . , Xj_i) ■ (yXi+i, . . . , Xj,-_i) ■ {Xj+i, . . . , Xn) 

= djdiAxi,x2,...,x„ {xi, . . . , Xj, . . . , Xn) ) 
they are well defined. Most of the time we will be interested in 

Axi,X2,...,Xn (yXi, X2, . . . , Xn) , 

which will be denoted simply by A (Xi, X2, . . . , X„). 

On the other hand, if the n-tuple {Xi, X2, . . . , X„} is fixed, define the polynomials A^ (x) by 

A^{xi : i = u{j) for some j) = ^x„(i),x„(2),...,x„(fe) {xu(i), Xu(2), ■ ■ ■,Xu(k)) ■ 
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This is a monic polynomial family. 

Example 7. From conditions (1211) and (l22t is follows that each A {Xi, . . . , X„) is a polynomial in 
Xi, . . . , X„ and their joint moments. Moreover, from conditions ^TH and (l23t it follows that this 
polynomial has degree 1 in each Xi. Some low-order polynomials are: 

A(Xi) =Xi-v9[Xi], 

A (Xi, X2) = X1X2 - V [Xi] X2 - <^ [X2] Xi - v5 [X1X2] + 2^ [Xi] [X2] , 

A (Xi, X2, X3) = X1X2X3 - [Xi] X2X3 - ^ [X2] X1X3 - (/^ [X3] X1X2 

- V [X1X2] X3 - [X2X3] Xi + [Xi] V9 [X2] X3 + [Xi] [X3] X2 

+ 2^ [X2] (/^ [X3] Xi - ^ [X1X2X3] + [X1X2] (/^ [X3] + [Xi] [X2X3] 

+ <^ [X1X3] ^ [X2] - 5^ [Xi] ^ [X2] ^ [X3] . 

Lemma 8. For n>l, the map A, 

(Xi, X2, . . . , X„) A (Xi, X2, . . . , X„) 

n-linear. 

Proof. For Xi = aX + we will show that 
(24) A«x+/3y,X2,...,x„ (aa; + X2, . . . , x„) 

= aAx,X2,...,X„ {X, X2, . . . , Xn) + PAy,X2,...,X^ {y, X2,..., Xn) ■ 

Then in particular, 

A{aX + (3Y,X2,...,X^)=aA (X, X2, . . . , X„) + (F, X2, . . . , X„) , 
and a similar proof holds for the other components. 

da:AaX+(SY,X2,...,X^ {<^X + (3y , X2, ■ ■ ■ , Xi, . . . , X„) = aAx2,...,X,, {X2, ■■■,Xn) 

= {aAx,x2,...,x^ (a;, X2, ■ ■ ■ , x„) + (3Ay,x2,...,x,, {y, 3^2, • • • , Xn)) ■ 

So condition (1^ implies that each monomial containing x in av4x,x2,...,x„ {x, X2, ■ ■ ■ , Xn) appears, 
with the same coefficient, in Aax+t3Y,X2,...,x„ (ctx + Py, X2, ■ ■ ■ , Xi, . . . , with x replaced by 
{ax + Py), and all monomials containing {ax + (3y) appear in this way. Application of dy pro- 
duces a similar statement, and implies in particular that the right-hand-side of (EU is a polynomial 
in {ax + I3y), and that monomials containing this term on both sides of (EH) coincide. By induction, 

diAaX+f3Y,X2,...,X^ {<^X + (3y,X2,...,Xi,..., Xn) 

= AaX+t3Y,...,X,.i {ax + f3y,. . Xi-i) ■ Ax,+ ,,...,X„ {Xi+l, ■■■,Xn) 
= {X,-.-, Xi_i) ■ Ax,+ ^,...,X„ {Xi+1, ■■■,Xn) 

+ /3^y,...,X._i (y, • • • , ■ {x^+l, ...,Xn) 

= di {aAx,x2,...,Xr, (a;, Xa, • . . , x„) + (3Ay.x2,...,x,, {y, X2, ■ ■ ■ , • 

This implies that all the terms on both sides of (l24l) which contain ax + Py, X2, . . . ,Xn coincide. 
Finally, 

^ [AaX+l3Y,X2,...,Xn (aX + /3F,X2,...,X„)] = 

= Lf [aAx,X2,...,X^ {X, X2, . . . , Xn) + PAy^X2,...,X^ {X, X2, . . . , Xn)] ■ 
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implies that the constant terms coincide as well. □ 

In the following power series, the x variables commute with the z variables, although neither the x 
nor the z variables commute among themselves. 

Lemma 9. Let X = (Xi, X2, . . . , X„) be a family of random variables, and R the corresponding 
free cumulant functional. Then in the expansion 

(1 - X ■ z + R(z))-^ = 1 + 5Z ^« W-^-' 

u 

Vc{l,...,n} ■K&IntiV) 
{TT,V)eNC(n) 

Proof. This is just formula (flTt for the coefficients of the inverse power series to 

n 

1 - Y^i^i - 



2=1 



□ 



-1 



Theorem 10. Let X = {Xi, X2, . . . , X„) be a family of random variables, and M, R, A be the 
corresponding moment and free cumulant functionals and free Appell polynomials. Then 

l + Y'^u (x) z^ = if(x, z) = (1 - X ■ z + R(z))-^ . 

u 

Proof. Using the substitution in Proposition^ 

i7(x,z) = (^l + i?(^u;i(l + M(w)),...,u;,(l + M(w))) - x • w(l + M(w))^ 

= (^l + M(w)-x-w(l + M(w))) ^ 
= (1 + M(w))"^(l - X- w)"\ 

Since 

(1 - X ■ w)"^ = ^ + X^Wit, 
u 

it follows that 

<f [H{X, z)] = (1 + M(w))-V [(1 - X ■ w)-^] = 1. 
So in the notation of Lemma |9l (p [i?t/(X)] = 0. 
Since 

(9^^(1 - X • z + R(z))-^ = (1 - X ■ z + R{z))-^Zi{l - X ■ z + R(z))-\ 
dxiBu = -Bu(i),...,MO-i)-BM(j+i),...,u(A;). 

j:u{j)=i 

Finally, it follows from Lemma |9l that in each monomial of -81,2,. ...n^ the variables appear in the 
increasing order of the indices. Therefore by the definition of free Appell polynomials, i?i,2,...,n = 

^l,2,...,n- 
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This is true for an arbitrary family of random variables, in particular for 

Using the explicit expression from Lemma |9l -Bu(x) is the free Appell polynomial Ai^2,...,k for Y 
applied to . . . , Xu(k)), which is exactly (for X). □ 

Corollary 11. The one-variable polynomials of Definition^are A^"\xi) = Aii i{xi). 

Proposition 12. Let u G W\ 

(a) The expansion of in terms of the free Appell polynomials is 

n fc+1 

k=0 Bc{l,2,...,n}, j=l 

(25) B={v{l)<v{2)<...<v(k)} 

= E (n^[^(-^)]+ E A{X^,.,s))llR[X^u:C)]), 

where v{0) = 0, f (A; + 1) = n + 1, 

(b) The recursion relation for the free Appell polynomials is 

n 

(26) A {X(^j,u)) = XjA {Xff) — ^ R[Xj, Xu(l), . . . , Xu(n~k)]A {Xu(n-k+l), ■ ■ ■ , Xu{n)) ■ 

k=0 

(c) The expansion of the free Appell polynomial is 

A{Xii)= ^ -^(uiV) ^ (-l)''''i?^[X(^^.yc)]. 

yc{l,...,n} vre/nUV^) 
i7r,V)eNC{n) 

Proof From the proof of the preceding theorem, 

(1 - X ■ w)-i = (1 + M(w)) + W ^«)' 

u 

where 

=Wi(l + M(w)). 

But 

(1 - X ■ w)"^ = ^ + E 

u 

So IS the coefficient of w-a in 

oo 

5^ 5^ A., (x) (1 + M(w))t/7,(i) (1 + M(w)) . . . (1 + M(w)) . 
fe=o v^nk 

The first line of equation (E51) follows. The second line, including the fact that B is outer, follows 
from the definition of the free cumulants. 

The recursion relation follows by expanding the identity 

(l-x-z + R(z))(l + ^A„^(x)2;„^) = 1. 

u 

The final expansion follows by combining the preceding theorem with Lemma|51 □ 
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3.4. Free binomial. Suppose tt g V{n) has the property that the collections {Xi : i E B}^^^ are 
freely independent. Then the basic relation between free independence and free cumulants iH^ 
says that 

(27) R(z) = ^R(zi : i G fi). 

BGtt 

More precisely, if R is the free cumulant generating function for {Xi, . . . , X„} (and so a function 
of zi, Zn), and for each subset i? C {1, n}, R^ is the free cumulant generating function for 
{Xi : i G B} (and so a function of {zi : i G -B}), then R = X^BeTr -^^s- 

Now let u G Afc „. It can be written uniquely as m = (vi,V2, . . . , Vk) with each Vi longest consecu- 
tive sequence from the same class of vr. 

Proposition 13. In the setting above, 

k 

A{X.a) = \[A{X,:). 
1=1 

The proof is based on the following lemma due to Franz Lehner. The lemma deals with bounded 
operators, but it applies equally well to formal power series. 

Lemma 14. [33 1 Let Si, . . . , Sn G B{H) be arbitrary operators and assume that the sum of 
alternating products 

oo 

/ + 5^ SiiSi^ . . . Si^ 

n=l 11^12^. 

(the sum over all products where neighboring factors are different) converges, then it equals 

-1 



^ N \ 

\ ?=1 / 



Proof of Proposition\n\ Let tt have classes B2, . . . , B^- From (l27t . 

TV 

(28) 1 - X ■ z + R(z) = 1 - XI ~ 



j=i ieBj 



Denote 



(l-X-z + R(z)) ' = 1 + ^ J2 ^i^u)z^ 

^==1 «e{l,...,n}* 



and for j = 1, . . . , A^, 



ieB J k=l t/e(Bj)ft 



24 M. ANSHELEVICH 

Using Lemma (141 with Sj = Hj — 1 and equation (l28l) . we get 



N 



1 + E E (^^1 - - 1) . . . - 1) = I 1 - - l)Hr^ 

n=l 



i=l 
N 



-1 



if. 



Equating coefficients of 2;^ on both sides of this equation gives the proposition. □ 
As a consequence, we get the "non-commutative binomial" formula: for X, Y freely independent, 
+ = A{X + Y,...,X + Y) 

n 

= J2Y1 [^^"^'^^ (X) ^("(2)) )^W3)) (X) . . . + (y ) (X) . . 

k=l MeAfe_„ 

This is the extension to the free Appell polynomials of the binomial expansion of (X + F)" for 
X, Y non-commuting variables. 

A different analog of the commutative binomial expansion (fT4b is the "co-multiplication" property 
of the free Appell polynomials: it easily follows from Definition |4l that for general x, y. 



n— 1 



x-y 



k=0 



where the x and 7/-dependent polynomials correspond to possibly different functionals. In particu- 
lar. 



n-l 



fc=0 

3.5. Martingale property. The martingale property of free Appell (and, more generally, Sheffer) 
polynomials for processes with freely independent increments was shown previously by Biane 1 1 5 1 . 
The following is an alternative proof for distributions all of whose moments are finite, which uses 
the binomial property above. 

Lemma 15. Let {^t} be a free convolution semigroup with all moments finite, {X{t)} the corre- 
sponding free Levy process, and jA^") j the corresponding free Appell polynomials. Then for each 
n, the process | A*^") (X(t)) } is a martingale with respect to the filtration induced by {X(t)}. 

We again note that A^'^\X{t)) = A^^'j^^{X{t)), and the polynomial ^^(t) depends on t, so the 



short-hand notation should be handled with care. 
Proof. 

E, [AW(X(t))] = E, (x{s) + {X{t) - X{s)) 
= E., 



^ J2 A("«)(X(s))AW2))(x(t)-X(s))... 



k=l ueAi, 



A(")(X(s)) 
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since E [A^''^ {X{t) - X{s))] = E [A(^) {X{s))] = for A; > 0, and using the definition of free 
independence and properties of conditional expectation. □ 

Theorem 16. Let Ui G W'^'^\ i = 1,2, . . . ,k, N = Yli=i -^l^)- ^ = ("^i; "^2, • • • , Uk). Each of 
the quantities 



(29) M[X^„X^„ X„^J = M[X^] 

(30) M[A {Xu,) , A {X^,) ,...,A (X„^J] 

(31) R[{X^,),{Xu,),...,{Xu,)] 

(32) R[A (X„^J ,A{X^,),...,A (X„^J] 



is equal to the sum of Rt^[X^ over partitions tt G NC{N) which are: 

(a) Equation (1291 ) ; arbitrary, 

(b) Equation ( 13 Oi l; non-homogeneous, 

(c) Equation (13 II) ; connected, 

(d) Equation (I32I) ; connected and non-homogeneous 

with respect to the interval partition 7!'s(i),s{2),...,siky Note that (BUI) Z5 f/ze linearization coefficient 
for the free Appell polynomials. 

There is a similar theorem for the usual Appell polynomials, obtained by replacing the lattice of 
non-crossing partitions with the lattice of all partitions. The following proof is also very similar to 
the one in [1251 . 

Proof. Equation (E^ is the basic relation between moments and free cumulants. Equation (BTT ) 
was proven in BO. i45 1 . We prove equation (1^ : the proof of ( BUt is similar, and also implied by 
equation (1^ below. In fact, we will prove a more general statement, that 

(33) R[A{X^,),...,A [X^;) , {X^^J, (X„^J] 

is equal to the sum of Rt,[X^ over all non-crossing connected partitions with no homogeneous 
classes in 113(1), s(2),...,s(j)- The proof will proceed by induction on s(l), . . . , s{j), starting with the 
statement for 

R[X^(i), Xy(^j), . . . , (X^^J] 

which is valid by equation (BTT) . Suppose ( B^ holds for all smaller j or, for the same j, for all 
shorter ui, . . . ,Uj. Substituting the recursion relation 

A (X,^) = X,^ - E E ^ i^i^f-B)) n Rl^i^r-c)] 

neNCisij)) Ben CStt, 

into equation (B^ breaks it up into the difference of two terms. The first term contains only (j — 1) 
Appell entries. It is the sum over all non-crossing connected partitions with no homogeneous classes 
in TXs{i),s(2),...,s{j-i)- The second term contains j Appell entries, the first (j — 1) of which have indices 
Ml, ... , Mj-i, and the last one is a proper sub-index of Uj. For each choice of a proper subset 

i?c|5^.(.) + l,...,^s(z)|, 

I 1=1 i=l ) 

the second term contains the sum over all non-crossing connected partitions such that the comple- 
ment of B in this subset consists of homogeneous classes, and tt has no homogeneous classes with 
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respect to tTs{i),s{2),...,b- In other words, the second term is the sum over all non-crossing connected 
partitions which have no homogeneous classes with respect to ns(i)^s{2),...,s{j-i) but have some ho- 
mogeneous classes in the j'th class of ns{i)^s{2),...,s{k)- Clearly the difference of these terms is the 
sum over all non-crossing connected partitions with no homogeneous classes in 'n's{i),s{'2),---,s(j)- '-' 

The following proposition is a generalization of (l30b . Its proof is similar. 

Proposition 17. Products of free Appell polynomials have an expansion 

(34) A{X^JA{X^,)...A{X^^)= Yl Yl RAX(u:V^^)]A{X^^.,v)). 

VC{1,...,N} TT&NCiV) 

non-homogeneous , 
in,V)eNCiN), 
V€ 0uter{TT,V) 

Here, "non-homogeneous" is with respect to the restriction 'ns{i),s{'i),-Ak) \ 

Proposition 18. Fix {Xi,X2, . . . ,Xn} with the moment generating function M, free cumulant 
generating function R, and free Appell polynomials {A^}. With respect to the joint distribution 
functional (/^x, 

(a) {A^} are pseudo-orthogonal iff{Xi, X2, . . . , Xn} form a semicircular family. 

(b) {Aft} are orthogonal iff{Xi, X2, . . . , Xn} form a free semicircular family. 

Proof. Suppose {A^} are pseudo-orthogonal. By formula (l30b . for k > 2, 

^ [A {Xu{i)) A (X„(2), . . . , Xu{k))] = R[Xu{i), Xu(2), • • • , Xu(k)] = 0. 

So all the joint cumulants of {Xi, X2, . . . , X„} of order greater than 2 are zero, and these random 
variables form a semicircular family. In particular, in formula (l30b only pair partitions make a non- 
zero contribution. If in addition {A^} are orthogonal, from this formula it follows in addition that 
R[Xi, Xj] = unless i = j. Therefore 

n 

R(z) = 5^(i?[Xi,X,,];22 ^ jiix,]z,). 

i=l 

So {Xi, X2, . . . , Xn} are freely independent, and form a free semicircular family. □ 

Remark 8. I thank Dima Shlyakhtenko for bringing to my attention the following observation. 
Given a family of self-adjoint elements {X^} whose free Appell polynomials are pseudo-orthogo- 
nal, there need not exist a linear change of variables B such that = Sj=i ^ijXj | is a family of 

self-adjoint elements whose free Appell polynomials are orthogonal. Indeed, one can always find 
such a linear transformation with R[Yi, Yj] = 5ij, but the Fj's need not be self-adjoint. As a result, 
if [Y*Yj] need not be 0. In other words, while the joint free cumulants of the 1^'s are 0, the joint 
free star-cumulants of the Fj's need not be 0, so the Fj's need not be freely independent. 

If (p has the trace property (p [ab] = ip [ba] , then one can always orthogonalize pseudo-orthogonal 
polynomials. 

3.6. Free Sheffer and Meixner families. Sheffer families are martingale polynomials for the cor- 
responding Levy processes. Based on this idea and the result of 1031 . in Q we proposed the 
definition of free Sheffer families, which are families of martingale polynomials with respect to 
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free Levy processes. Specifically, free Sheffer families are families of monic polynomials whose 
ordinary generating function H{x, t, z) = J2'^=o Pn{x, t)z^ has the form 

1 1 1 



l + f(z)t-Uiz)x l + fiz)t 1 



U{z) 



■X 



l+f{z)f 

We also described all the free Meixner families, that is, free Sheffer families consisting of orthogo- 
nal polynomials. They are given by the recursion relations 

Pi(x, t) = X, 
xPi{x, t) = P2{x, t) + ax + t, 

xPn{x,t) = Pn+l{x,t) + aPn{x,t) + (t + 6)P„_i(x,t) 

for n > 2. For a = b = 0, the polynomials are the Chebyshev polynomials of the second kind. For 
general a G M, 6 G M+, such polynomials have been considered by a number of authors; see the 
discussion in pages 26-28 of [ 12], and also ll24ll . 

Definition 6. Let R be a free cumulant generating function, and U(z) an n-tuple of non-commuta- 
tive power series such that Ui{z) = Zi+ higher-order terms. Multivariate free Sheffer polynomials 
are defined via their generating function 

if(x, z) = (1 - X ■ U(z) + R(U(z)))-' = l + Yl Pu{^)^u- 

u 

Define a linear operator Dj (left partial derivative) on non-commutative power series via 

Theorem 19. Suppose that for a family of multivariate free Sheffer polynomials, Pq and Pi are 
orthogonal for all i and \u\ > 2. In particular, this is the case if the polynomials are pseudo- 
orthogonal. Suppose also that the covariance matrix R[Xj, Xi] is non- degenerate. Then 

n n 

U(z) = F<-i> i?[Xi, X,]z„ . . . , ^ P[X„, X,]z^ . 

i=l 1=1 

Here Fi{z) = (DjR)(z) — ip [Xi], and for an n-tuple of power series F, G = F<^^^ is the inverse 
of F under composition, 

F,(G(z)) = Zi. 

In the classical case, the corresponding condition defines precisely the natural exponential families. 
The proof below is inspired by the one in f^H^. 

Proof Using the substitution Ui{z) = Wi{l + M(w)) leads to R(U(z)) = M(w). Then 
i/(x,z) = (^1-x- w(l + M(w)) +M(w)) ^ = (l + M(w))"^(l-x-w)"\ 

It follows that [H(X., z)] = 1. Suppose the conditions of the theorem are satisfied. By definition, 

Pj = Xj-^ [X,j\, so P; = Pj. Then 



^[p;z,H{y.,7)\=^ 



n 

E 

i=l 



if [PjPi] ZjZi. 
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So 



and 



n 

R[Xi, X,]zjZ, = ^ [PjZjH{X, z)] = if [XjZjH{X, z)] - ip [Xj] 



i=l 



ip [XjZjH{X, z)] = J2 Mh^i + ^ [^j] 
j=l 

On the other hand, using the substitution above, 

XjZjH{x, z) = XjZj(l - X ■ U(z) + R(U(z)))^"^ 



Zj. 



and so 



Since 



^XjZj{l-x-w{l + M{w)) + M{w)) ^ 

= zj(l + M(w))-S.(l - X • w)-^ = zj(l + M(w))-^ (xj + J2 

u 

[XjZjH{X, z)] = zj{l + M{w))-' (^M[Xj] + M[Xj,Xa]ws^ . 

it 

DjM{w) = Dj (J2 M[Xa]w^^ = M[Xj] + ^ M[Xj,Xa]wa, 



it follows that 

<p [XjZjH{X, z)] = zj{l + M(w))-^D,(M(w)). 
Since Ui(z) = ^^(1 + M(w)) and R(U(z)) = M(w), 

(l + M(w))-iD^,(M(w)) = (Z^,^R)(U(z)). 



We conclude that 
and 

Thus 

Since 



ip[X,z,H{X,z)] = Zj{D,K){V{z)) 

n 



1=1 



F,(U(z)) = J]i?[X,-,X,]^,. 



Fj{z) — ^ R[Xj, Xi\zi + higher-order terms 



i=l 



and the covariance matrix is non-degenerate, this series has an inverse under composition. 
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3.7. Free Kailath-Segall polynomials. We only list three results in this section; all the other re- 
sults from Section have free analogs, which will be proven in greater generality in Section 
Free Kailath-Segall polynomials are the g-Kailath-Segall polynomials from that section for q = 0. 

First, free Kailath-Segall polynomials have explicit expansions 

TTeInt{n) ScSing(7r) BeS C£S'= 

Second, the product of free Kailath Segall polynomials can be expanded as 

k 

n W {fu,{l), fu,i2), fu.isii))) 
i=l 

TT&NC{N) Sen, BeS'= 

'rA7r,(i),,(2),...,,(,)=oSmg(5-)=0 

(the proof is similar to the one in Proposition [l]). Third, free Appell polynomials can be expressed 
in terms of the free Kailath-Segall polynomials. 

Proposition 20. For xi = x, 

n 

^^"n:^) = E E 

k=l ueAk.n 

Proof. The proof is similar to Proposition|5l We show that 

n k 
k=l nSAfc^^ i=l 

by induction, with the induction step being to verify that 

n+l r- k k k k - 

A(-^)(x).i=$: E -ill 

k=l ugAj. i=l i=2 i=2 i=l 

n j k 

-E^"+wE E H^-ii)- ° 

j=0 k=l veAkj i=l 

4. g-INTERPOLATION 

We saw in the preceding two sections that both the Appell and the Kailath-Segall polynomials 
are related to probability theory (they are martingale polynomials), and that they have analogs 
performing the same functions in free probability theory. On the other hand, there is also a very 
clear relationship between these classes and the third one of orthogonal polynomials. In the light of 
this, it is interesting to look at their g-deformations. For orthogonal polynomials, such deformations 
are given by the members of the Askey scheme of basic hypergeometric orthogonal polynomials 
ll29ll . On the other hand, a possible g-deformed probability theory was initiated by Bozejko and 
Speicher [16|. We show that the Kailath-Segall polynomials can be extended to this context, with 
the same relation to the deformed probability theory and orthogonal polynomials as in the q = 
0, 1 cases. As a consequence, we obtain combinatorial formulas for some basic hypergeometric 
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orthogonal polynomials, and Wick product formulas for Fock representations. In contrast, the q- 
extension of the Appell polynomials is not satisfactory. We propose a definition in the one- variable 
case which fits well with the orthogonal polynomials picture (the analogs of the Meixner families 
are the Al-Salam and Chihara polynomials). However, this definition does not fit well with the 
g-deformed Levy processes (see the Appendix), and an adequate definition in the multivariate case 
is missing. 

4.1. Notation. Fix q E (—1, 1). A few standard pieces of g-notation are: [0]g = 0, 



forn > 1, [0],! = 1, H,! = UtM,iorn> 1, 
Define the g-cumulants for a single measure via the relation 

Define the cumulant generating function for this section to be 

oo ^ 



n=l 



More generally, for a n-tuple of random variables, define the joint g-cumulants via 

TTeNC{n) -BGtt 

4.2. g-Appell polynomials. There is a number of different possible definitions for g- Appell poly- 
nomials. For example, two g-deformations of the relation dxAn = nAn~i, leading to such defini- 
tions, were considered in |IT1|4|. We prefer, instead, to use an interpolation between the recursion 
relations in the classical and the free cases. 

Definition 7. The g- Appell polynomials are defined via the recursion relation 

An+l{x) = xAn{x) - X] K j rn+l-kAk{x), 

k=0 ^ ^ 9 

where {r^} are some g-cumulant sequence. 

Proposition 21. The generating function ofq-Appell polynomials has the form 

oo ^ oo ^ 

g MJ^"^"^^"" = ^^^''^ = Jll- (l- q)xzq^ + R{zq^) -Rizq^^^y 
Proof Denote Dg^^z" = [nlqz""^, so that 

z — qz 

is the standard g-derivative operator. By definition, 

1 1 _ " _ 1 1 

A _ „ ^ A \^ ^ ^ A n+1 
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Divide by [n + 1]^ to get 
Apply Dq^^ to get 

1 1 " 1 1 



So 

It follows that 



Hix.z) — r ^ , , =— — -H(x,qz). 

^ ' ^ l-{l-q)xz + R{z)-R{qz) ^ ^ 



□ 



Example 9. Since orthogonal polynomials satisfy three-term recursion relations, the only orthogo- 
nal polynomials among the g-Appell are those with = for k > 2. Thus Il{z) = az + bz'^. By 
adding a constant to x and re-scaling the polynomials, we may take = 62kt, R(^) = j^^'^' 
which case 

Up to further re-scaling, this is the generating function for the continuous (Rogers) g-Hermite poly- 
nomials. 

Example 10. Consider the polynomial family defined by the recursion 

xPn = Pn+1 + {t- [n]q)Pn + [n]qXPn-l. 

It has the generating function satisfying 

H{x,t,z) = — — — ^- ^ — H{x,t,qz). 

i + (i-?)t^^- {i-q)zx 

So this is a family of g-Appell polynomials for R(2;) = t [nf^^' ^ ~ ^^9- As g — >^ 0, 

H converges to 

1 

1 + j^t-zx' 

1—z 

the generating function for the free Appell polynomials for the free Poisson process. On the other 
hand, for g ^ 1, H converges to 

e^"(l - zf, 

the generating function for the Appell polynomials for the Gamma process. This can be considered 
as an extension of a well-known property that for t = 1, a square of a normal random variable 
has a Gamma distribution, while a square of a semicircular random variable has a free Poisson 
distribution. 

Proposition 22. Let {A^} be the q-Appell polynomials for the q-cumulant generating function R. 
Then 

n-l 



k=0 ^'^^1 
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where {A^} are defined via their generating function 



?i=0 

Proof. By Proposition l2T 



An\\ ^ l-{l-q)xz + Ii{z)-n{qzy 



d,H{x, z) n 1 _ (1 _ q),:zqk + _ r(^^/c+i) 

oo ^ 
" 2 ^(^' ^) 1 _ (1 _ + R(^gfe) _ R(^gfc+1) " ^^^^^^^ 

OO oo oo 

n=0 ^ fc^o [iJ<?- 

CO ^ CO ^ ^ 

= E m E iry^'^^^T^TTT^'^' 

n=0 1=0 ^ ^ 

oo oo ^ 

EE m „,,„, A.M-4?M."^". □ 



n=0 i=0 L J? L ' J? 

Example 11. If {P„(x)} are the continuous g-Hermite polynomials, then 



n— 1 r 

\n 



k=0 

where 



d^Pn{x) = -^—-Pk{x,t)Qn-k-l{x), 



E'5"' 



XI 



„^o 1 - (1 - q)xz + (1 - g)^^ 

So Qn{x)z"' = f/„(-\/l — gx)(v/l — g-z)", where f/„ are the Chebyshev polynomials of the second 
kind. Therefore 

d.Pnix) = J2r^^^AV^)'''Pk{x)U^-k-i{VT^qx). 

4.3. Al-Salam-Chihara polynomials. In this section we consider a g-deformation of the Meixner 
families. One such deformation has been considered in 1 27 1. It is based on the exponential function 
for the operator calculus for the Askey-Wilson operator [26|. Previously, Al-Salam showed that 
under this approach, the unique orthogonal Appell family also consists of (multiples of) Rogers 
g-Hermite polynomials Under the more elementary approach of Q, the unique orthogonal 
Appell family are the Al-Salam-Carlitz polynomials. 

However, we have an extra requirement to put on our families. In addition to the correct limiting 
behavior as g — > 1, we also require a correct limiting behavior as g — > 0. Specifically, in this limit 
we should obtain the free Meixner families. 

Define orthogonal polynomials by the recursion relation 

(35) xPn = Pn+1 + a[n]gP„, + [n]g{t + b[n - l]g)P„_i 
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for a E M., b E M+ (cf. [7. Remark 6]). Then the free Meixner families are these families of 
polynomials for g = 0, while the classical Meixner families are these families for g = 1. 

Lemma 23. Up to re-scaling, the polynomials defined by recursion (1351) are the Al-Salam-Chihara 
polynomials. 

Proof. Define the generating function H{x, t, z) = Yl^=o l^i^nix, t)z^. Then 

xz{H{x, t, z)-l) = - Pii^^ t)^ - H{x, t, qz) + Pi{x, t)qz] 

OjZ 

+ ^—^[H{x,t,z) - H{x,t,qz)] 
bz^ 

+ tz^H{x, t, z) + -[H{x, t, z) - H{x, t, qz)]. 

Since Pi{x, t) = x, 

, / 1 az r, bz"^ \ N / 1 CLZ bz'" 
H{x, t, z) xz + h + = H{x, t, qz) h h 



1 — g 1 — g 



Therefore 



\ (1 + az + bz'^) , 

H{x,t,z) = — , , 2^ I n ^771 TH{x,t,qz) 

[I + az + bz'^j + (1 — q)[tz^ — xz) 

(l + az + bz^) 

H (x, t, qz) 



1 + [a - (1 - q)x\z + [6 + (1 - q)t\z'^ 

H{x, t, qz). 



1 



Affine transformations w = Jb + (1 — q)tz, y = g)^ ° ^ bring H to the standard form of the 
generating function for the Al-Salam-Chihara polynomials 11291 3.8] 



rj! / ,/N 0''w){l-b'w) , , 

H{y, w;a,b) = — ■ ^H{y, qw;a,b), 

^ 2yw + 



where a' , b' are the roots of the quadratic polynomial z"^ H — , " z + b. □ 

y'b+(l-q)t 

In particular, the Hermite case a = b = corresponds to the Rogers (continuous) g-Hermite 
polynomials, while the Charlier case 6 = corresponds to what are usually called the continuous 
big g-Hermite polynomials. Note also that the g-Krawtchouk polynomials defined in iHTIl are, up to 
a shift, of this form, with a = 1 — 2p, b = —p{l — p),2t = Np{l — p), except that in this case, as 
expected, b is negative. 

Al-Salam-Chihara polynomials were defined in fl] as all polynomials other than the Meixner fam- 
ilies characterized by a certain convolution property. Their measure of orthogonality was found 
explicitly in [11 \. The interpretation of these polynomials as g-analogs of the Meixner families was 
explicitly conjectured in [5|, and proved in Q. We found the following proof independently, and 
in our particular case it is also somewhat simpler. For a more interesting characterization using 
stochastic processes, see lITHI . 
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Theorem 24. Suppose the generating function H{x, z) = J2'^=o '\^.^n{x)z'^ of the monic orthog- 
onal polynomials defined by the recursion relation 

has the form 

Hix, z) = F(z) |T , , = — - — ———H(x, qz), 

lJ-Q^-i^-(l)Uiq''z)x F{qz)l-{l-q)U{z)x ^ ^' 

where F,U are formal power series with F{z) = 1+ higher-order terms, U{z) = z+ higher-order 
terms. Then the polynomials are a re-scaling of the Al-Salam-Chihara polynomials, with 

and 

Pn = [n]q{Pi + C2[n - l]g). 

Proof By assumption, 

(1 - (1 - q)Uiz)x)H{x, z) = 4r^H{x, qz). 

Therefore 

(36) .H(., .) = [m., .) - f^^H(., . 

Define the lowering operator D on M.[x] by D(P„) = [n]gPn-i and extend linearly. Then D{H) = 
zH and equation (l36t implies 

1 / Fiz) 
D{xH){x, z) = I zH{x,z) — — - — -qzH{x,qz) 



(37-, :i - q)Uiz) V ' ' ' F{qz) 

Expand jjhs^r^ into the formal power series X^^o ^it^ cq = 1. Equation dTTl says 



(38) D{xPn){x) = [n],xP„_i(x) + V Cfc^^^g"-'^P„_fc(x). 

Applying the lowering operator to the recursion relation gives 

(39) D{xPn)ix) = [n + l]gP„(x) + a„[n]gP„_i(a;) + /3„[n - l],P„_2(a;). 
The recursion relation forn — 1 is 

xPn-l{x) = Pn{x) + an-lPn-lix) + /3„_iP„_2 (x) . 

Multiplying it by [n]g gives 

(40) [n]gxPn-i{x) = [n]gPa{x) + [n]ga„_iP„_i(x) + [n]gPn-iPn-2ix). 



APPELL POLYNOMIALS AND THEIR RELATIVES 35 

Combining Equations (l38t . (l39t . and (l40b . we obtain 

- [n + l]gPn - a„[n]gP„_i - /?„[n - l]qP„_2 = 0. 

Collecting coefficients of P„_i gives 

(41) + ci[n]gg""^ - a„[n]g = 0; 
collecting coefficients of Pn-2 gives 

(42) + C2[n]Jn - - /5„[r2 - 1], = 0, 
and collecting coefficients of Pn-k for A; > 2 gives 

Cfc = for /c > 2. 

Equation (HTT) gives 



while equation (B^ gives 



n— 1 

a„ - a„_i = cig 



(3n—l n—2 



We conclude that a„ = + Ci[n]g and = [n]g(/3i + C2[n — □ 
Thus it appears reasonable to define the g-Sheffer polynomials via their generating function 

oo ^ oo ^ 

^^^^ S M?^"^'"' ^ ]?o 1 - (1 - g)xf/(;.g'=) + K{U{zq^)) - R(f/(zg'=+i)) " 

4.4. g-Kailath-Segall polynomials. The origin of these polynomials is in the g-Levy processes 
defined in [6, 9J. 

Definition 8. Let be a complex star- algebra without identity, and (■) a star- linear functional on it. 
Let A be the complex unital star- algebra generated by non-commuting symbols {X{f) : / G ^o''} 
(and 1) subject to the linearity relations 

X{af + Pg) = aXU)+PX{g). 

Equivalently, A is the tensor algebra of ^o- The star-operation on it is determined by the require- 
ment that all X{f), f e Aq"" are self-adjoint. For such /j, define the g-Kailath-Segall polynomials 

byW(/)=X(/)-(/)and 

n 

W(/,/i,/2,...,/n)=X(/)W(/i,/2,...,/n)-$^g^-M//.)w(/l,...,/,,...,/. 

(44) 

-Y^q'-'W (//,,...,/„...,/„) -(/)W(/i,/2,...,/„). 

1=1 

This map has a C-linear extension, so that each W is really a multi-linear map from Aq to A. 



36 M. ANSHELEVICH 



n+1 '^-^ ~ " n-k 



In the particular case = Co[x] (polynomials without constant term), we may denote Xi = X(x*). 
The functional can be taken to be the g-cumulant functional of a measure ji, {x'') = ri{fi). Then 

are multivariate polynomials in {xi : i E N}. 

Proposition 25. The Kailath-Segall formula (fT5t takes the following form. 

W{f, /, ._J) = J2(-lf^JS—Xif'^-'')W{f, f,...J)-{f) W{f, /,...,/> 



[n-k],V 

The proposition can be proven directly by induction, or deduced from Theorem l26l b). 

Many of the following formulas appear in [21] in the g-Gaussian case = Co(x) and (a;^) = 5|u|,2- 

Theorem 26. The following expansions hold. 

(a) A monomial in X{f) 's can be expanded in terms of the q-Kailath-Segall polynomials: 

X(/i)X(/2) . . . X(/0 = ^'^^''^^ n W (/c : C G ^) . 

(b) For a permutation a G Sym(n), write its standard cycle decomposition as 

a = (m(1,1),...,m(1,Si)) ... {u{k,l),...,u{k,Sk)). 

Here 1) = min^ j)}, and the cycles of a are ordered according to the order of their 
smallest elements, 1) < u{2, !)<...< m(/c, 1). Let Sing((T) be the one-element cycles 
of a, and let s{a) the number of inversions of the permutation 



Fia) 



1 ... n 

m(1,1) ... u{k,Sk) 



(F is almost, but not quite, the fundamental transformation of Foata |23|j. Finally, for a 
subset S C Sing(cr), denote 

-^(5,o-)(/l, f2, ■ ■ ■ , fn) = Y\. X{fu{i,l)fu{i,2) ■ ■ ■ fu{i,Si))- 

i:{u{i,l),...,u{i,Si))GS'^ 

Note that in fact, X(^s,cr) depends only on {fi : {i} ^ S}. With this notation, the q-Kailath- 
Segall polynomials are 

w(/i,/2,...,/„)= E (-ir^^^^'^^+i^ig^^'^^ n (/.)^(5,.)(/i,/2,...,/„). 

(T€Sym(?i) 5cSing(a-) i:{i}£S 

Proof. Part (a) was proven in the appendix of |9| using the Fock space representation; one can also 
use the defining recursion relation. For part (b), by definition, 

n 

W(/,/i,/2,...,/n) =X(/)W(/i,/2,...,/n)- 5^ g*-M//.)W (/!,...,/„...,/„) 

i=l 

n 

- J^g-^W (//„...,/„...,/„) -(/)W(/i,/2,...,/n). 



i=l 
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Using the defining recursion and induction on n, this is 

TeSym(n) 5cSing(T) jeS 

n 

- E (//^) E E n (z^) ^(5>-)(/i' •••,/.,•••,/«) 

»=1 TGSym({l,i,n}) 5cSing(T) j€5 



n 



»=1 TeSym({o,l,i,n}) ■'^'^ 

n 



^=1 TeSym(|0,l,i,n» ■5cSing(r) jeS 

X X(^s,T){ffi, fl, ■ ■ ■ , fi, ■ ■ ■ , fn) 

-(f) H E (-ir^'^^^^^'"5^^^^n(/^)^(5.-)(/i'- ••'/") 

TeSym(n) 5cSing(T) je5 

Since in the fourth term, e Sing(r), the second and the fourth terms cancel. So we obtain 

E E (-Ij^'^+^^-^^^^^^+^^+l^lg^^^) n if^) ^if)^iS,r){fl, ■■■Jn) 

TeSym(n) ScSing(T) jeS 
n 

+ E E E n (Z^) /i, ■■■,/.■■■ , /n) 

'=1 TeSym({0,l,i,n}) '^^^™|('^) 

TeSym(n) 5cSing(T) jeS 

The three terms in the preceding equation correspond to pairs (5*, cr), a e Sym ({0, 1, . . . , n}) such 
that 

e Sing (a) ,0 ^ S (first term), 

^ Sing(a), (t(0) = i (second term), 

Oe S C Sing (a) (third term). 

It remains to match up the powers of q. Suppose that 

a = {u{l, 1), . . . , u{l, si)) . . . {u{k, 1), . . . , u{k, Sk)), 

s{a) = i(F(a)) for 

^^^)= («(1,1) u{k,s,))- 

In the first and the third terms. 



Fia) 



1 ... n 

u{2,l) ... u{k,Sk)J ' 
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and i{F (a)) = i(F(r)) for 

F(r) 

In the second term, 

and i{F{a)) = z(F(r)) + i - 1 for 

F(t) - 



1 ... n 

u(2, 1) ... u{k,Sk] 

1... n 

i ... u{k,Sk] 



1 ... n 

3) or m(2, 1) ... u{k,Sk) 



Define a unital linear functional on ^ by 

^[W(/i,/2,...,/n)]=0 

for all n > 0, /i, . . . , /„ G A 

Corollary 27. The functional ip is given by 

(45) ^ [X(/i)X(/2) . . . X(/„)] = 5^ g-W n (f^) ■ 

■KeV(n) Be-K 

That is, with respect to ip, the q-cumulants of such an n-tuple are 

R[X{f,),X{f,),...,X{f^)] = {fj2...fn). 

Corollary 28. The linearization coefficient for the q-Kailath-Segall polynomials is 



□ 



n W (/n,{l), /n,{2), • • • , fuMi))) 



,i=l 



E """" n ■ 



'rA7rg(i)_j,(2)_....s(fe) = 
Sing(7r)=0 



In other words, it is the sum of partitioned q-cumulants over all inhomogeneous partitions with no 
singletons, with weights q^^^'^\ 

4.5. Fock space realization. The g-Kailath-Segall polynomials form a monic polynomial family 
in infinitely many variables. So we can construct a Fock space for them as in Section ITSl this space 
will be infinite-dimensional. Instead, we will construct the Fock space directly from the multi-linear 
maps W. As a vector space, it will be the space of all polynomials in elements of Aq, modulo the 
linearity relations. Equivalently, it is -^o The induced inner product is determined by 

I k n \ n 

(46) ((g)Aw,(8)/.0))='5n. E 9^^'^^n(/-0)A'^0))>- 

\ i=l j=l I (7eSym(n) j=l 

If (•) is a faithful state on Aq, for q E (—1, 1) this inner product is positive definite ifT^ . Hence 
is a state. 

For the classical case q = I, this inner product is only positive semi-definite even if (•) is positive 
definite. The quotient by the kernel of (p gives precisely the symmetric Fock space. 
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From the defining recursion relation, the action of the operator X(/) is 

® /2 ® • • • ® /„) = / ® /l ® /2 ® • • • ® /„ 

n 

i=l 

n 

1=1 

+ (/) /l ® /2 ® • • • ® U 

Thus it is a sum of a creation operator 

a*if) : /l ® /2 ® • • • ® /n / ® /l ® /2 ® • • • ® /n, 

an annihilation operator 

n 

a{f) : fi® f2® ...® fn^Yl fl® ...® fi® ...® fn, 

1=1 

a preservation operator 

n 

P{f) ■■ fl® f2® ...® fn^J2 I'^^ff^ ® . . . ® /^ ® . . . ® /n, 

1=1 

and a scalar operator 

/l ® /2 ® . . . ® /„ ^ (/) /l ® /2 ® • • • ® /n- 

Now consider the case when (■) is positive semi-definite but not faithful. For simplicity, we will 
also assume that it has the trace property, (fg) = (gf) for all f,gE Aq. Note that this does not 
imply that (/? is a trace. Denote by 

/(.) = {/ G A : (/7) = 0} 

the kernel of (■). Denote by Hq the Hilbert space obtained by completing the vector space A/I(^.) 
with respect to the norm induced by ||/|| = (/*/), with the induced inner product. Suppose /(.) is 
in fact an ideal, so that 

(47) V/ e I^.^Wg e Ao, {rg*gf) = 0. 

This is the case, for example, when Aq is a C* -algebra. Note that 

\{f9)\<V{f*f) (9*9), 

and so (^AqI^.)) = 0. It follows that the Fock representation of A factors through to the representa- 
tion on 0^^o nf'' with the inner product induced from (l46b . 

Note that Ao/I(-) is an algebra. In the examples below, we will observe the following situations. 
Let CI © ^0 be the standard unital extension of the non-unital algebra ^o- 

(a) If the functional (■) can be extended to CI © ^ in a positive way, then | (/) |^ < (/*/) (1), 
and so (/( )) = 0. Thus the action of X(f) in the Fock representation depends only on the 
class of / in Ao/I{.)- So we may replace A^ by A//(.) throughout. 
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(b) If the functional (■) has no such positive extension, it is natural to take X{f) to be the sum 
of only three operators, X(f) = a*(f) + a(/) +p{f). The formulas for the g-Kailath-Segall 
polynomials need to be modified accordingly; this approach was used in |9|. Under this 
construction, again the action of X(f) in the Fock representation depends only on the class 
of / in Ao/I(.), and so we may replace by Ao/I(.) throughout. 

(c) A particular case of part (b) is when 

(48) AoAo C /(.). 

In this case in the representation of Ao/I{.), X{f) = a*{f) + a(f). 

Example 12. Throughout this example, = Co[x]. 

(a) Let be a positive measure on M all of whose moments are finite, and (x*) = mj(z/) for 
i > 1. Then we are in the context of part (a) of the preceding alternative, and so may replace 
^0 by Ao/I{-) throughout. This is the case of compound g-Poisson distributions. 

More specifically, let (x*) = tfor i > 1, so that u = t5i. Then 

/(.) = {P : P(l) = 0} 

is an ideal. This is the g-Poisson case. In Aq /!{.), = x for all i. So the g-Kailath-Segall 
polynomials in this case can be considered as single- variable polynomials. Specifically, 
they are the continuous big g-Hermite polynomials; the recursion relation for their centered 
version is a particular case of equation (l35t . Theorem l26l and Corollary |28] provide combi- 
natorial identities for these polynomials. 

(b) Let be a positive measure on M with all moments finite, and (x) = 0, (x*) = mj_2('^)- 
Then (■) is positive on Co[x], but in general has no positive extension to all of C[x]. So 
we are in the context of part (b) of the preceding alternative. As a result, we may replace 
^0 by Aq/Ii^.) throughout, as long as the operators X(x*) without the scalar part are used. 
This is the g-Kolmogorov case (see which includes all the centered g-infinitely divisible 
measures (all of whose moments are finite). 

(c) Let ^0 = Co[x], (x*) = t5i^2 = t'^i-2{^o)- Note that this functional is positive on Co[x], but 
not on C[x]. Then 

/(.) = {Pe Co[x] : P'(0) = 0} 

is an ideal, and condition is satisfied. This is the g-Gaussian case. In Aq/Ii^.), x* = for 
i > 2. \i follows that the g-Kailath-Segall polynomials in this case can also be considered 
as single- variable polynomials. Specifically, they are the continuous g-Hermite polynomials 
of Example |3 Theorem EH and Corollary |2H] provide combinatorial identities for these 
polynomials as well. 

4.6. Differences from the g = 0, 1 cases. Binomial families of polynomials are Sheffer families 
for t = 0. In particular, for any Appell polynomials the binomial family is always {x"}. 

In Q, we investigated g-analogs of binomial families, which have generating functions 

CO ^ 

(49) exp^([/(z)x) = n 1 WTT^- 

We showed (Proposition 18) that the theory of such binomial families for various values of g is 
exactly parallel to the theory considered by Rota et al. for g = 1. In particular, the lowering 
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operator (in the sense of Theorem l24h for such a sequence is U^^{Dq). In fact, these results were 
surely known before. 

However, the following is one difference between the classical and free Sheffer (or even Appell) 
families. For a classical (single variable) Sheffer family {P„(a;,t)} with generating functions 

(50) U-\d^)Pn{x,t) = Pn-i{x,t). 

However, for a corresponding free Sheffer family with generating function , the corre- 

sponding property is 

n-l 

U~\d,)Pn{x,t) = ^Pfe(x,t)P„„fc_i(x,t). 
k=0 

In other words, the lowering operator for the classical family is U^^{dx) independently of t. In the 
free case, for t = the lowering operator is U^^{Do), but it changes for t > 0. 

For g 7^ 0, 1, the binomial families corresponding to the g-Sheffer sequences of this paper, with 
generating functions 

oo ^ 

n 1 _ U{q^z)x' 

k=0 ' 

are different from the binomial families ( l49b . For binomial families of this paper, the lowering 
operator need not be a function of Dq (which is the lowering operator for the family {a;"}). For 
example, the binomial family for the continuous big g-Hermite polynomials is 

n-l 
k=0 

and the lowering operator for this family is 

f{x)-f{l + qx) 
X — (1 + qx) 

This operator does not commute with Dq. 

Another property which still holds in the free case is that any free Sheffer family can be expressed 
as a linear combination of the corresponding Appell family |I71 Lemma 1]. This property does not 
hold in the q case: see the appendix. 

Appendix A. g- Appell polynomials are not martingale polynomials eor g-LEVY 

PROCESSES 

All the calculations below are performed with Maple 7. 

A. 1 . Relation between g- Appell and g-Sheffer polynomials. Let {P„ (x, t) } be the centered con- 
tinuous big g-Hermite polynomials, and t)} the g- Appell polynomials for ri = 0, = t for 
k > 1. Note that {Pn} are g-Sheffer for this cumulant sequence. Then 

P4(x, t) - (Ai{x, t) - (3 + 2g + q^)A3{x, t) + (3 + 4g + 3q^ + g^) ^(x, t) 

- (l + 2g + 2g2 + g3)Ai(x,t)' 
= t[xg(g - 1)]. 



42 



M. ANSHELEVICH 



It follows that g-Sheffer polynomials need not be a t-independent linear combination of the corre- 
sponding g-Appell polynomials. In particular, since {P„(x,t)} are known to be martingale poly- 
nomials for the g-Poisson process LSJ, it follows that the g-Appell polynomials, in this case, are 
not. 

A. 2. Generic case. Take to be the polynomial algebra generated by symbols {x(t) : t E M+}, 
and the functional on it be given by 

(x(ti)x(t2) • • = (mintj)rfc. 

i 

Denote X{t) = X{x(t)) as in the case (b) of Section |431 Then {X{t)} is a g-Levy process in the 
sense of [i6J. For each fixed t, define Wii{t) using the recursion relation (l44b involving only the 
for that specific t. The conditional expectations onto the subalgebras generated by the 
g-Levy process {X(t)} are determined by 

E. [W^{t)] = W^{s). 

In particular, for each u the Kailath-Segall family Wn^t) is a martingale. In a number of situations, 
these processes also have single-variable martingale polynomials families. That is, for each n there 
is a family of polynomials {P{x, t)}^^^ oo) °f degree n'mx such that 

(51) E,[P{X{t),t)] = P{X{s),s). 

This is the case for g = 1, g = 0, if the process is a g-Brownian motion, and if the process is a 
g-Poisson process. In all of these cases, from the existence of such martingale polynomials one can 
deduce the Markov property for the corresponding process [il7n 8ll. 

We show that generically, there is no degree 5 polynomial which is a martingale for such a g-Levy 
process. This is a strong indication that general g-Levy processes do not have the Markov property. 
For a fixed t, explicit Maple calculations (see the author's web page) show that there is a monic 
degree 5 polynomial P^i^x), whose coefficients depend only on t, g, and the free cumulants of the 
process, such that 

P,{X{t)) + q\l - q)t{rsW2{t) - r^W^it)) 

= time-independent linear combination of the Kailath-Segall polynomials. 

Since a linear combination of Kailath-Segall polynomials is a martingale, the expression above is a 
martingale. Suppose there is a degree 5 martingale polynomial in x. Further explicit calculations 
show that there definitely are martingale polynomials of degrees 4 and less. It follows from equa- 
tion ( BTT ) that the leading coefficients of such polynomials have to be independent of t, so we may 
assume them to be monic. By subtracting from P5 the degree 5 martingale polynomial, we obtain a 
martingale 

P4(X(t)) + q\l - q)t{r3W2it) - rsW^aW), 

where the polynomial P4 has degree at most 4. P4(X(t)) ■ 1 contains a term of degree degP4, 
while {r3W2(t) — r2W3{t)) ■ 1 = {r^x{tY — r2x(t)^), so if degP^ > 1, P4 also has to have a 
time-independent leading coefficient. Proceeding in this fashion, we obtain a martingale of the 
form 

Y{t) = a{t, g, r)X{t) + b{t, q, r) + q''{l - g)t(r3W^2(t) - r2Ws{t)) 

Since 

Y{t) ■ 1 = a{t, g, r)x(t) + b{t, g, r) + q\l - q)t{r^x{tf - r2x{tf), 
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b is independent of t, and may be taken to be zero. First assume a = 0. This can occur only in the 
following cases. 

(a) q = 1- This is the classical Levy process case. 

(b) g = 0. This is the free Levy process case. 

(c) r^W2{t) = r2W3{t). Suppose that condition ( BTl) holds. From the "creation" part of the 
representation of the operator X(t) in Section l43l it follows that r3x{ty — r2x(tY G /(.), 
and in particular that 

{x{tf{r^x{tf - r^xitf)) = t{r^rk+2 - rar-fc+a) = 

for > 0. This is the case only when = a^^"^ for k > 2 (we may take ri = 0, r2 = 1). 
Then Wk{t) = a''~^X(t). These are g-Poisson processes with step a, and the degenerate 
case a = gives the g-Brownian motion. 

For general a, and g 7^ 0, 1, 

E, [Y{t)] ■ 1 = ait, q, r)x(s) + g'(l - qW^xis)^ - r2x(s)=^), 
and the martingale condition implies 

{a{t,q,r) - a(s, g, r))x(s) +g^(l - q){t - s)(r3x(s)^ - r2x{s)^) = 0. 
As above, this implies 

(a(t, g, r) - a{s, g, r))rk+i + g^(l - q){t - s){r^rk+2 - '^a'^fc+s) = 

for /c > 0. By fixing s and varying t, it follows that a(t, g, r) is linear in t, and may be taken to be 
a(g, r)t. With r2 = 1, rs = a, this means that 

a 

rk+3 = otrk+2 + -irr^ ^^k+l■ 

r(i - q) 

For general a, a, this says that {r^} is a sum of two geometric sequences, corresponding to the 
g-Levy measure being supported at two points. 

Corollary 29. The Markov processes of ifTSi are not q-Levy processes. 

Proof. Proposition 3.3 of 1.1 8J shows that the processes of that paper have martingale polynomials. 

□ 
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